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We present an all-sky formalism for the Cosmic Microwave Background (CMB) bispec- 
trum induced by the primordial non-Gaussianities not only in scalar but also in vector and 
tensor fluctuations. We find that the bispectrum can be formed in an explicitly rationally 
invariant way by taking into account the angular and polarization dependences of the vector 
and tensor modes. To demonstrate this and present how to use our formalism, we consider 
a specific example of the correlation between two scalars and a graviton as the source of 
non-Gaussianity. As a result, we show that the CMB reduced bispectrum of the intensity 
anisotropies is evaluated as a function of the multipole and the coupling constant between 
two scalars and a graviton denoted by gtss', \beee\ ~ x 8 x 10~^*|(?tss|. By estimating 
the signal-to-noise ratio, we find that the constraint as \gtss\ < 6 will be expected from the 
PLANCK experiment. 

§1. Introduction 

Bispectrum (three-point correlation functions) of the Cosmic Microwave Back- 
ground (CMB) temperature anisotropies has been attracting attention over the years 
as a powerful observational tool to investigate the primordial non-Gaussianity.'^''^ 
As is well known, primordial curvature perturbations whose statistics deviate from 
the pure Gaussian ones can produce the nonzero bispectrum of the CMB temperature 
anisotropies. The size of the primordial non-Gaussianity has often been parametrized 
by a so-called nonlinearity parameter /nl- Depending on the shape of the bispec- 
trum, this nonlinearity parameter can be sorted into three types of /nl, called local, 
equilateral and orthogonal types. Current observational limits on these /nls are 
given by -10 < f^l-"^ < 74 for the local type, -214 < /^'f " < 266 for the equi- 
lateral type and —410 < f^l^°^ < 6 for the orthogonal type (95% C.L.).^' These 
observational constraints are still consistent with the Gaussian primordial curvature 
perturbations that are expected to be generated from the standard single slow-roll 
inflation model. However, one can also consider that there may be some cue for hunt- 
ing the non-Gaussianity because the central values of some /nls have deviated from 
zero. Hence, if future experiments would confirm that the statistics of the primordial 
curvature perturbations deviate from the Gaussian ones, then the standard single 
slow-roll inflation model can be excluded as a dominant mechanism of generating 
primordial fluctuations. Thus, the primordial non-Gaussianity can be considered as 
a new probe of the mechanism of generating primordial curvature perturbations. 

As discussed above, previous studies have focused largely on non-Gaussianities 
in scalar-mode perturbations. However, we can also consider the sources of non- 
Gaussianities in vector and tensor perturbations, for example, the nonlinear cou- 
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plings between gravitons and scalars during inflation,'^' nonlinearities of the Sachs- 
Wolfe effect P'EJS 'EI cosmic stringsP'liS and primordial magnetic &elds^^^^^^^ 
Hence, in order to constrain the size of non-Gaussianities and understand the na- 
ture of these sources, one should include the contribution of vector and tensor non- 
Gaussianities to the CMB bispectrum. 

In our previous work,1Sl' we presented the bispectrum formulae of the CMB 
temperature and polarization anisotropies sourced from non-Gaussianity not only in 
scalar but also in vector and tensor fluctuations. It is expected that the nonlinearities 
will also induce the tensor and vector mode perturbations, and the modes may 
generate more characteristic features in the CMB angular spectra than in the scalar 
one. In Ref. [TB|) , we found that the bispectrum formulae for vector and tensor 
modes in all sky analysis formally take complicated forms compared with the scalar 
mode case owing to the dependence of the photon transfer functions on the azimuthal 
angle between the wave vector of photon fluctuation k and the unit vector specifying 
the line of sight direction n. However, by using flat sky approximation, we have 
simplified the equations of bispectra of the CMB anisotropies to solve the above 
difficulty because no azimuthal dependence arises in this limit. 

This paper is an extension of our previous work. We present a general formalism 
of the CMB bispectrum induced from the primordial vector and tensor fluctuations 
in the all-sky analysis. We newly consider the angular dependences in the polariza- 
tion vector and tensor bases, which have been neglected in our previous work. To 
demonstrate how to calculate the CMB bispectrum by making use of our formula, 
we show a calculation of the CMB bispectrum induced from the primordial non- 
Gaussianity generated through the interaction between two scalars and a graviton 
(tensor) during inflation, which has been originally discussed by Maldacena.® 

This paper is organized as follows. In the next section, we present a formulation 
of the CMB bispectrum in the all-sky approach, which is an extension of our previ- 
ous work."^ In ^ we show the calculation of the CMB bispectrum induced from 
the primordial non-Gaussianity generated through the two scalars and a graviton 
(tensor) correlator during inflation. We define a coupling parameter characterizing 
the strength of such interaction as gtss and evaluate an observational limit on gtss 
by calculating the signal-to-noise ratio. In the final section, we give a summary and 
conclusion of this paper. 

§2. Formulation of the CMB bispectra for scalar, vector and tensor 

modes 

In this section, we derive general formulae of the CMB bispectrum of temper- 
ature and polarization fluctuations induced by the primordial non-Gaussianity in 
scalar, vector or tensor-mode perturbations in the all-sky analysis. 

At first, we introduce an expression of CMB fluctuation. In the all-sky analysis, 
CMB fluctuations of intensity or polarization field are expanded with the spin-0 or 
spin-2 spherical harmonics, respectively .^^"^'^ Then, the coefficients of CMB 
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fluctuations, cafled a^m, are described as 



M-if j ^5;[sgn(A)]^+-_Al^4(fc)e(')(fc)ri5(fe) , (2-1) 



where the index Z denotes the mode of perturbations: Z = S (scalar), = V (vector) 
or = T (tensor) and its hehcity is expressed by A; A = for {Z = S), = ±1 
for [Z = y) or = ±2 for [Z = T), X discriminates between intensity and two 
polarization (electric and magnetic) modes, respectively, as X = I,E,B and x is 

determined by it: x = for X = I,E or = 1 for X = B, ^^^^ is the initial perturbation 

(z) 

decomposed on each helicity state and 7J ^ is the time-integrated transfer function 



in each sector (calculated in, for example, Refs. [T8|) .[2T ]) .[22 ]) )r)| 
Next, we expand S,^^^ with spin-(— A) spherical harmonics as 



CWik)^Y.i'2ikUYUk) , (2-2) 



£m 

and eliminate the angular dependence in Eq. (|2-1|) by performing fc-integral: 



A 

Here, we use the orthogonality relation of spin-A spherical harmonics a^^SJ,!!!) 

d'^hxYl,^,{h)xYira{n) = 6i^i>6m,m' ■ (2-4) 
The initial bispectrum in vector or tensor-mode perturbations will be expressed 



This definition, which includes the angular dependence on k in the polarization vector 
or tensor, is more general than Eq. (3) of Ref. [T8|) . We will see in the discussion 
in ^that the initial bispectrum from inflation indeed takes the above form for the 
tensor case. On the other hand, to calculate the CMB bispectrum using Eq. (|2-3p . 
the bispectrum of is needed. If the primordial bispectrum satisfles the rotational 
invar iance, we can set it as 

(nd:,i(^o)-(2vr)^-^i^^53'^ X) ■ (2-6) 

Then, using Eq. ()2-2p . the conversion equation between F^^^^As g^j^j ^^^^^h 
rived as 



mim2m3 



*■* Here, we set = 1. 
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(2-7) 



From Eqs. (\'2-'d\i . (|2-6p and the orthogonality of Wigner-3j symbols as Eq. (|A-5p . the 
CMB angle- averaged bispectrum, which is defined a^'l^ 

^(ZiZaZg) = f ^1 ^2 ^3 

XiXiX-i/i/i/'i ^ \ rrii 1712 ms 




can be written as 



Tj{ZiZ2Zi) 



n=l ^^^^ A„ 



x{27TfTl^l^j^-'ik^,k2,k3) . (2-9) 

Thus, when one computes the CMB bispectrum, only the alternative initial bispec- 
trum J^^^i^fi'^'^ is necessary in each case. 

§3. CMB bispectrum induced by the primordial non-Gaussianity in the 

two scalars and a graviton correlator 

In this section, we demonstrate how to calculate ^£^£^}^^ the CMB bis- 
pectrum by considering the contribution of two scalars and a graviton correlator. SI 
Furthermore, we evaluate an observational limit on the primordial non-Gaussianity 
of the graviton sector by calculating the signal-to-noise ratio. 

3.1. Two scalars and a graviton interaction during inflation 

We consider a general single-field inflation model with Einstein-Hilbelt actiorPSJ ; 



S= d^x 



(3-1) 



where g is the determinant of the metric, R is the Ricci scalar, M^^ = l/(87rG), (/> is a 
scalar field, and X = —g'^'^d^(j)du(t)/'2,. Using the background equations, the slow-roll 
parameter and the sound speed for perturbations are given by 

_ H _ Xp^x 2 _ P,x , . 

' - mMl, ' - 2Xp,xx +P,x ' ^ 

where H is the Hubble parameter, the dot means a derivative with respect to the 
physical time t and p^x denotes partial derivative of p with respect to X. We write 
a metric by ADM formalism 

ds^ = -N'^dt^ + a^e^-^idx" + N^dt^dx^ + N''dt) , (3-3) 

where and A^" are respectively the lapse function and shift vector, ^ab is a trans- 
verse and traceless tensor as 'jaa = da^ab = 0, and e'''"'' = 5ab + lab + ladcbl'^ H ■ 
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On the flat hypersurface, the gauge-invariant curvature perturbation C, is related to 
the first-order fluctuation of the scalar field as ^ = —Hip/ (p. Following the con- 
version equations (|B-20p and (|B-27p . we decompose C, and ^ab into the helicity states 
as 

= C(fc) , i^^^\k) = \e^^^\k)^ab{k) . (3-4) 

(±2) 

Here, e^j^ is a transverse and traceless polarization tensor explained in Appendix 
iBl The interaction parts of this action have been derived by Maldacena^ up to the 
third-order terms. In particular, we will focus on an interaction between two scalars 
and a graviton. This is because the correlation between a small wave number of the 
tensor mode and large wave numbers of the scalar modes will remain despite the 
tensor mode decays after the mode reenters the cosmic horizon. We find a leading 
term of the two scalars and a graviton interaction in the action coming from the 
matter part of the Lagrangian through X as 

^Isrd-order ^ a~'^^^labda^dbip , (3-5) 

therefore, the interaction part is given by 

5'mt 3 j (fxagtsslabdaCdbQ. (3-6) 

Here, we introduce a coupling constant gtss- Prom the definition of Cilab and the 
slow-roll parameter, gtss = £• For a general consideration, let us deal with gtss as a 
free parameter. In this sense, constraining on this parameter may offer a probe of 
the nature of infiation and gravity in the early universe. The primordial bispectrum 
is then computed using in-in formalism in the next subsection. 

3.2. Calculation of the initial bispectrum 

In the same manner as discussed in Ref . , we calculate the primordial bispec- 
trum generated from two scalars and a graviton in the lowest order of the slow-roll 
parameter: 

Xeif ^(fcl)fc2a46 , (3-7) 

jfT T T ^ , , kik2 + + ksh kik2k3 

I{ki,k2,k3) = -kt^ r \ (3-8) 

where kt = ki + k2 + k^^ and * means that it is evaluated at the time of horizon 
crossing, i.e., a*i7* = k. Here, we keep the angular and polarization dependences, 
^ab (^i)^2a^36) which have sometimes been omitted in the literature for simplic- 
ity.'^''^''^ We show, however, that expanding this term with spin-weighted spher- 
ical harmonics enables us to formulate the rotational-invariant bispectrum in an ex- 
plicit way. The statistically isotropic power spectra of ^'^'^'^ and are respectively 
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given by 

^ {27rfPs{mk - k') 

k^Psjk) _ Hi 

^ {2^f^^5(k - k')5x,x' (for A = ±2) , 
= ^ = 8cs.e.As = -^5 , (3-9) 

where r is the tensor-to-scalar ratio and As is the amphtude of primordial curvature 
perturbations. Note that the power spectra satisfy the scale invariance because we 
consider them in the lowest order of the slow-roll parameter. Using these equations, 
we parametrize the initial bispectrum in this case from Eqs. ()3-7p and (|2-5p as 

F±2»»(S,„S,2,S;3) = /"''"(ti,t2,»:3)eif '(fci)t2At , (3-10) 

n,.kM.'-!^^^^^t^^. (3.11) 

12 3 ^ 1 

Note that f^'^^^^ seems not to depend on the tensor-to-scalar ratio. In Fig. [1] , we 
show the shape of I/ki. From this, we confirm that the initial bispectrum /C^'^'^) 
()3-lip dominates in the squeezed limit as ki <^k2 — k^ like the local-type bispectrum 
of scalar modes. 

In the squeezed limit, the ratio of /(^'^'^) to the scalar-scalar-scalar counterpart 
j{SSS) _ ^fj^i^Pg(^ki)Ps{k2), which has been considered frequently, reads 

f(^^10gts^I_ktk2 20gtss I , . 

/(^^^) 3/nl h kl 3/nl h- ^ ■ ^ 

In the standard slow-roll inflation model, this ratio becomes 0(1) and does not 
depend on the tensor-to-scalar ratio because gtss and /nl are proportional to the 
slow-roll parameter e, and I/kt has a nearly flat shape. The average of amplitude 
is evaluated as I/kt ^ —0.6537. Therefore, it manifests the comparable importance 
of the higher order correlations of tensor modes to the scalar ones in the standard 
inflation scenario. 

3.3. Formulation of the CMB bispectrum 

For this case, by substituting Eq. ()3-10p into Eq. ()2-7p . the initial bispectrum is 
given by 

3 

rr 

x^2Y,\^^{k^)Y,l^^{k2)Y,l^^{ks) 

x/(^^^)(A:i,A:2,A:3)eif)(fei)4a465 (^J^fei^ • (3-13) 
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Fig. 1. (color online) Shape of 7/fci. For the symmetric property and the triangle condition, we 
limit the plot range as k\ < k2 < fca and |fci — fc2| < fcs < fci + ^2. 



We derive this simpler form as the following procedure. 

At first, we express all parts containing the angular dependence with the spin 
spherical harmonics: 



3b 



Mmarrif, 



M rUa mi) 



(3-14) 



\i=i 



llY.{-l)^^/^jL^{hy)Yl,j^ik, 

i=l LiMi 

Li L2 L3 



rO 

'^^L^L2L3 1 j^^ j^^ 



where we used the relations listed in Appendices [A] and |B] and 



JS1S2S3 



(2/i + l)(2/2 + l)(2/3 + l) / h h h 
47r V S2 S3 



(3-15) 



(3-16) 



Secondly, using Eq. (|A-8p . we replace all the integrals of spin spherical harmonics 
with the Wigner symbols: 



j d'fci T2Yl,rnSkr)Yl^M,{ki)±2YiM{ki) = /JfJ 
<fk2 Yl^^,^{k2)Yl.^M2ik2)Y{^Ak2) = Ill'i 



d'ks Y;^^^{ks)Yl^M,{k3)Y{„.,{ks) = 



2{ h Li 2 

mi Ml M 

i2 L2 1 

m2 M2 ma 

h ^3 1 

m-s M3 mt 



(3-17) 
,(3-18) 
.(3-19) 
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Thirdly, using the summation formula of five Wigner- 3j symbols as Eq. (lA^ . we 
sum up the Wigner-3j symbols with respect to azimuthal quantum numbers in the 
above equations and express with the Wigner-9j symbol as 



M1M2M3 

Mmami, 



Li 
Ml 



L2 
M2 



Li 
Ms 



2 

M 



1 



1 

rub 



mi 




1 


)( 


is 


rria 






h 


h 




Li 


L2 




2 


1 


1 



L3 
M3 



1 

nib 



(3-20) 



After these treatments, performing the summation over mi, 7712 and 7713 like Eq. (|A-5p . 
we can obtain the final form as 

(87r)3/2 
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-f(^''\ki,k2,ks) 



ErO 
-'LiL2L3-'€iLi2 

L1L2L3 



t±20t2 tO p 



rO 



«1 

Li 
2 



«2 
1 



4 
1 



y^dy 



(3-21) 



Note that the absence of the summation over nii,m2 and in this equation means 
that the tensor-scalar-scalar bispectrum maintains the rotational invariance. As 
described above, this consequence is derived from the angular dependence in the 
polarization tensor. Also in vector modes, if their power spectra obey the statis- 
tical isotropy like Eq. ()3-9p . one can obtain the rotational invariant bispectrum by 
considering the angular dependence in the polarization vector as Eq. (jB-12p . 

Then, substituting the expression ()3-2ip into Eq. (|2-9p . we can calculate the 
CMB bispectrum induced from the nonlinear coupling between two scalars and a 
graviton. The CMB angle-averaged bispectrum is derived as 
JTSS) 



B 



(87r)3/2 



3 



L1+L2+L3 



7-0 7-20-2 tO 7-0 
^LiL2L3^£iLi2^i2L2l-^e3L3l 



Li 
2 



«2 
L2 
1 



L3 
1 



y'^dy 



i^x^ Jo 



f(^^'\ki,k2,ks) 



where we use the summation over Ai = ±2 as 

lAi+i'i 7-A1O-A 



9 7-20-2 

^-'^iLi2 



=±2 



(for xi + Li + £1 = even) , 
(for xi + Li + £1 = odd) . 



(3-22) 



(3-23) 
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Considering the selection rules of the Wigner symbols explained in Appendix [Aj we 
see that the bispectrum (j3-22p has nonzero value under the conditions: 

Li = |"'^^'"'^^ (forX,=/,i.) L3 = K3±1|, 

\\h±l\ (forXi =S) ' ' ' ' ^ ' ' ' 

3 

\Li - L2\ < L3 < Li + L2 , '^Li = even , 

i=l 

\h -i2\<h<h+h, yi. = = ' (3-24) 

' ^1 - ^ - i \odd (for Xi = B) . ^ ^ 

In Figs. [2] and [3l we describe the reduced CMB bispectra of intensity mode 
sourced from two scalars and a graviton coupling: 

ATSS) ASTS) ASST) 



- veihk) y^iii/ie2i3^^iii,eie2h^ iii/ihhj ' 1-3^0; 
and primordial curvature perturbations: 

"111/112^3 ~ VtMal ^IIIMf-2h ■ 

For the numerical computation, we modify the Boltzmann Code for Anisotropics 
in the Microwave Background (CAMB).'^''^ In the calculation of the Wigner-3j 
and 9j symbols, we use the Common Mathematical Library SLATEC^^' and the 
summation formula of three Wigner-6j symbols (|A-2ip . As the radiation transfer 
functions of scalar and tensor modes, namely, Tx\. and T^\,, we use the Ist-order 
formulae as discussed in Refs. I22p and I27p . From the behavior of each line shown 
in Fig. [3] at small £3 that the reduced CMB bispectrum is roughly proportional to 
we can confirm that the tensor-scalar-scalar bispectrum has a nearly squeezed- 
type configuration corresponding to the shape of the initial bispectrum as discussed 
above. From Fig. [2l by comparing the green dashed line with the red solid line 
roughly estimated as 

\ATSS) ASTS) ASST) i .-4 ^ o ^ 10-18 U I 97^ 

\"iii,m + "iii,m + "iii,m\ ~ * x 8 x iU \gtss\ , (3-27j 

we find that \gtss\ ~ 5 is comparable to /^l^' = 5 corresponding to the upper bound 
expected from the PLANCK experiment. In the next subsection, we check the 
validity of the above evaluation by computation of the signal-to-noise ratio assuming 
the zero-noise data. 

3.4. Estimation of the signal-to-noise ratio 

Here, we compute the signal-to-noise ratio by comparing the intensity bispec- 
trum of Eq. (j3-22p with the zero-noise (ideal) data and examine the bound on the 
absolute value of gtss- The formulation of (the square of) the signal-to-noise ratio 
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Fig. 2. (color online) Absolute values of the CMB reduced bispectra of temperature fluctuation for 
ii = £2 = £3- The lines correspond to the spectra generated from tensor-scalar-scalar correlation 
given by Eq. ()3-25|l with gtss = 5 (red solid line) and the primordial non-Gaussianity in the scalar 
curvature perturbations with /nl^' = 5 (green dashed line) . The other cosmological parameters 
are fixed to the mean values limited from WMAP-7yr data reported in Ref. [S]). 



(S/N) is reported in Refs. [T]) and [2]). In our case, it can be expressed as 



2 I n(TSS) d(STS) o(-55T) 



n) ~ ^ ^2 ' (3-28) 

' 2<£x<l2<l3<l hhh 

where cr^^^j^a denotes the variance of the bispectrum. Assuming the weakly non- 
Gaussianity, the variance can be estimated a^2Sli>l3S} 

^fi^2^3 ^ Ci^Ci^Ci^Ai^i^i^ , (3-29) 

where A^^i^n^ takes 1,6 or 2 for ^1 7^ ^2 7^ ^3,^1 = ^2 = ^3) or the case that two 
^'s are the same, respectively. denotes that the CMB angular power spectrum 
included the noise spectrum, which is neglected in our case. 

In Fig. m the numerical result of Eq. ()3-28p is presented. We find that (S/N) 
is a monotonically increasing function roughly proportional to i for i < 2000. It is 
compared with the order estimation of Eq. ()3-28p as Ref. [2]) 



2 



5^ ^j{2£y 



N ) V 24 V 47r 

-5 I 







,ATSS) ASTS) ASST), 

\"iiim "I" "mm "iiiim 



^lxbAxlQ-^\gtss\ . (3-30) 

2 



Here, we use Eq. (|3-27|) and the approximations as ^ ~ ^^/24, £^ 







0.36 X ^, and ^Cn ~ 6 x 10 We confirm that this is consistent with Fig. HI 
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Fig. 3. (color online) Absolute values of the CMB reduced bispectra of temperature fluctuation 
generated from tensor-scalar-scalar correlation given by Eq. (|3-25p (TSS) and the primordial 
non-Gaussianity in the scalar curvature perturbations (SSS) as a function of £3 with £1 and £2 
flxed to some values as indicated. The parameters are flxed to the same values deflned in Fig. [J] 



T 1 — I I I I I I 1 1 1 — I I I I I I 1 1 1 — I I I I I I 




0.0001 r 
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_i I ' I I I ' I I I 



70 



700 
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Fig. 4. (color online) Signal-to-noise ratio normalized by gtss as a function of the maximum value 
between £i,£2 and £3, namely, £. Each parameter is flxed to the same values deflned in Fig. [J] 



which justifies our numerical calculation in some sense. This figure shows that from 
the WMAP and PLANCK experimental datap^'l^ which are roughly noise- free at 
£ < 500 and 1000, respectively, expected {S/N)/gtss values are 0.072 and 0.16. 
Hence, to obtain (S/N) > 1, we need \gtss\ > 14 and 6. The latter value is consistent 
with a naive estimate \gtss\ ^ 5, which was discussed at the end of the previous 
section. 
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§4. Summary and discussion 

In this paper, we present a full-sky formalism of the CMB bispectrum sourced 
from the primordial non-Gaussianity not only in the scalar but also in the vector 
and tensor perturbations. As an extension of the previous formalism discussed in 
Ref. [TH]) . the new formalism contains the contribution of the polarization vector 
and tensor in the initial bispectrum. In Ref. [T8|l . we have shown that in the all- 
sky analysis, the CMB bispectrum of vector or tensor mode cannot be formed as a 
simple angle-averaged bispectrum in the same way as that of scalar mode. This is 
because the angular integrals over the wave number vectors have complexities for 
the non-orthogonality of spin spherical harmonics whose spin values differ from each 
other if one neglects the angular dependence of the polarization vector or tensor. In 
this study, however, we find that this difficulty vanishes if we maintain the angular 
dependence in the initial bispectrum. 

To present how to use our formalism, we compute the CMB bispectrum induced 
by the nonlinear mode-coupling between the two scalars and a graviton.® The 
typical value of the reduced bispectrum in temperature fluctuations is calculated as 
a function of the coupling constant between scalars and gravitons gtss- l^z/ffi^ + 

^\^iim + ^m£«l ~ X 8 X 10^^^\gtss\- Through the computation of the signal-to- 
noise ratio, we expect a constraint as \gtss\ < 14 from WMAP and \gtss\ < 6 from 
PLANCK. Although we do not include the effect of the polarization modes in the 
estimation of gtss in this study, they will provide more beneficial information of the 
nonlinear nature of the early universe. 

Our formalism will be applicable to the other sources of vector or tensor non- 
Gaussianity, such as, the cosmic stringd^''^® or the primordial magnetic fields.'^ 
Actually, in the specific case of vector-vector-vector correlation, we have already 
presented the rotationally invariant bispectrum from primordial magnetic fields.'^ 
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Appendix A 

Useful Properties of the Wigner Symbols 

Here, we briefly review the useful properties of the Wigner-3j, 6j and 9j symbols. 
The following discussions are based on Refs . [32 ]) . [331) . [Mil .[35 ]) . [361) . 

A.l. Wigner-?, j symbol 

In quantum mechanics, considering the coupling of two angular momenta as 

h = h + h, (A-1) 

the scalar product of eigenstates between the right-handed term and the left-handed 
one, namely, a Clebsch-Gordan coefficient, is related to the Wigner-3j symbol: 

h h ^3 \ _ (-l)^i-^2+"'3 {hmrl2m2\{hl2)hm^) 

This symbol vanishes unless the selection rules are satisfied as follows: 

\mi\ < li , \m2\ < h , l^isl <h, mi + m2 = ms , 

\h — h\ ^ h ^ h + h (the triangle condition) , li + I2 + h ^ ■ (A-3) 
Symmetries of the Wigner-3j symbol are given by 



h h 

mi m^ m2 



y mi m2 m^ J \ ^2 w-i f^s J 

(odd permutation of columns) 

h h h \ _ f h h h 
m2 m^, mi J y 777-3 "^1 ^2 

(even permutation of columns) 

Y —iT^i —m2 —ms 
(sign inversion of 7711,7712,7713) . (A' 

The Wigner-3j symbols satisfy the orthogonality as 

(2^3 + l)X] ( mi m2 mj, ) f mL mL m-, ) = '^"^i,™i'^™2,m^ ' 
\ h h I3 

•mim2 



mi 7772 ?7l3 
3 



h 


I2 


h 


m'l 


m'r. 


7773 


h 


12 




mi 


m2 





(2/3 + 1) > ^ ^ '1 '3 ^ , s ^, . (A-5) 



For a special case that Ya^^ k = even and ttii = 7772 = 7773 = 0, there is an analytical 
expression as 

h h h 




^3 (EL I) V(-^i +h + k)Wih -I2 + h)Wih + I2- k)l 
_l)E,=i — .V 1 — (A-6) 
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This vanishes for Ylf=i ~ odd. The Wigner-3j symbol is related to the spin- 
weighted spherical harmonics as 

2 

n.ii...(n)= s.yu(^)ia-r'' ( X X X ) ' ^^-'^ 

which leads to the "extended" Gaunt integral including spin dependence: 
Here ri''2S3 ^ / (2ii+i)(2/2+i)(2/3+i) f h h h 

hhh V 47r \^ _5j^ 52 S3 

A. 2. Wigner-Qj symbol 

Considering two other ways in the coupling of three angular momenta as 

^5=^1+^2+^4 (A-9) 

= h + h (A-10) 
= /i+/6, (A-11) 

the Wigner-6j symbol is defined using a Clebsch-Gordan coefficient between each 
eigenstate of Z5 corresponding to Eqs. (|A-10p and (|ATip as 

h h h \ _ (-i)'^+'^+'"+'M(^i^2)/3;/4;/5m5Mi;(/2 /4)/6;^5rri5) ^^,^2) 



Z4 Z5 Iq 



V(2/3 + l)(2/6 + l) 



This is expressed with the summation of three Wigner-3j symbols: 
(l)T.Lih-m, ( h h h \ 



X 



(A-13) 



7714 msmg 

I2 h \ I h h h 

h h ^3 \ f ^1 h h 
mi m2 w-s y \ U k h 

hence, the triangle conditions are given by 

|/i - h\ <h<h+ h, \k -k\<k<k + k , 

\h - k\ <k<h+ k, \k - h\<k<k + h ■ (A-14) 

The Wigner-6j symbol obeys 24 symmetries such as 

h h h\^{h h ^3 If h h h 
k h h ] \ h k h \ \ h h k 

k h ^3 1 _ f h h h 
h h h \ \ k h h 
(exchange of two corresponding elements between rows). (A- 15) 



(permutation of columns) 
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Geometrically, the Wigner-6j symbol is expressed using the tetrahedron composed 
of four triangles obeying Eq. ()A-14p . It is known that the Wigner-6j symbol is 
suppressed by the square root of the volume of the tetrahedron at high multipoles. 

A. 3. Wigner-9j symbol 

Considering two other ways in the coupling of four angular momenta as 

k = li + h + h + h (A-16) 
= l3 + l6 (A-17) 
= l7 + ls, (A-18) 

where Z3 = Zi + l2,le = Z4 + '5,^7 = Zi + ^4,^8 = I2 + I5, the Wigner 9j symbol 
expresses a Clebsch-Gordan coefficient between each eigenstate of Iq corresponding 
to Eqs. ([A471) and (|A48l) as 




V(2/3 + 1) (2/6 + l)(2/7 + l) (2/8 + 1) 
This is expressed with the summation of five Wigner-3j symbols: 



(A.19) 



E 

m4m5m6 
rriTrngmg 



k 



k 

7715 



^6 
TUq 



h 



k 
ms 




^3 



(A-20) 



and that of three Wigner-6j symbols: 



/i 
h 
h 



h 
h 




l)'"(2x + l) 



h 

X 



h 
k 



k 

X 



k 



h 

X 



k 
h 



19 
h 



(A-21) 



hence, the triangle conditions are given by 



1^1-/2! <l3<h + l2 , \k 
\h-k\ <h<h + k , Ih 



h\ ^ le ^ k + h 
h\ <h<h + h 



h\ <h<h + h , 

h\ <l9<k + k •(A-22) 



The Wigner-9j symbol obeys 72 symmetries: 
= (_l)Eti'« 





(_l)Eti^« 
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(odd permutation of rows or columns) 

h h h 
h h h 
h h h 
(even permutation of rows or columns) 





(reflection of the symbols) . (A- 23) 

Appendix B 

Polarization Vector and Tensor 



We summarize the relations and properties of a divergenceless polarization vector 
ei^^'' and a transverse and traceless polarization tensor e^^^^^^ . 

The polarization vector with respect to a unit vector n is expressed using two 
unit vectors 6 and <j) perpendicular to n as 

e^^^\n) = ^[Uh)±iMn)] ■ (B-1) 

This satisfies the relations: 

n'^6i±i)(n) = 0, 

eW(n)6f )(n) = (for A, A' = ±1) . (B-2) 

By defining a rotational matrix, which transforms a unit vector parallel to the z-axis, 
namely z, to n, as 

(cos On cos — sin sin On cos 4>n 
cos 9n sin (j)n cos 4>n SlU 9^ s'm (f)n I , (B-3) 
— sin On cos On 

we specify and cf) as 

e{h) = S{h)x , 0(n) = S{n)y , (B-4) 

where x and y are unit vectors parallel to x- and y-axes. By using Eq. (|BTp . the 
polarization tensor is constructed as 

eif)(n) = x/26i±i)(n)e[±^)(n). (6-5) 

To utilize the polarization vector and tensor in the calculation of this paper, we 
need to expand Eqs. (|B-ip and (|B-5p with spin spherical harmonics. An arbitrary 
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unit vector is expanded with the spin-0 spherical harmonics as 



-m{6m,i + (Jm,-l) 



(B-6) 



Here, note that the repeat of the index imphes the summation. The scalar product 
of a™ is calculated as 



/ 47r / 47r 

m m /T/'"<T/'" * X 

"a "a ~ 3 ^ ^ ^m—m' i '^a ~ ^ 'Jm,m' ■ 

Through the substitution of Eq. ()B-4p into Eq. ()B-6p , 6 is expanded as 

m mm' 
^ m 

Here, we use the properties of the Wigner D-matrix g^^^^^^Mi 



(B-7) 



47r 



1/2 



i)^_,y;„(n) 



21 + 1 

In the same manner, cf) is also calculated as 

(i)a{n) = + -l) X] «!r3"i^lm(«. 

^ m 

hence, the explicit form of Eq. ()B-ip is calculated as 

4^'^(n) = T5]a^±iYim(n). 



(B-8) 

(B-9) 
(B-IO) 

(B-11) 
(B-12) 



Substituting this into Eq. (|B-5p and using the relations of Appendix |X] and 



= 27^) the polarization tensor can also be expressed as 



This obeys the relations: 



eif)(n)=n,eir(n) = 0, 

"aft 



(B-13) 



elb' (^)elb ' (^) = 2<^A,- A' (for A, A' = ±2) 



(±2), 



(B.14) 
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Using the projection operators as 

O^e*'^-^ = k-^Vae'^-"^ = ikae'^-"^ , (B-15) 

Om,^k■. ^ (^fc-2v^v, + ^) e^*'- = (^-Kk + ^) e^^- , (B.16) 

Q^±i)gifc.a, ^ _i4±i)(fc)e*fc-^ , (B-17) 

Q(±i)gifc.., = fc^^VaOj^^^e*'^-^ = ^ae[^^^(fc)e^*=-^ , (B-18) 

the arbitrary scalar, vector and tensor are decomposed into the heUcity states as 

r]{k)=r]^^\k), (B-20) 

w„(/e) = a;W(fc)O(0) + ^ t^W(fc)OW , (B-21) 

A=±l 

= x^'Hk)0^^ + X^^\k)0i^ + X^'Hk)0l^^ . (B-22) 

A=±l A=±2 

Then, using Eq. ()B-2p and ()B-14p . we can find the inverse formulae as 

J'^\k) = -ikaOJaik) , (B-23) 

J^'\k) = ie^^'\k)ua{k) , (B-24) 



(0) (fc) = 3 ^ ^^^^^^ ^ ^g_25) 

X^^'Hk) = kae'i;'\k)xab{k) , (B-26) 
X(^')(fc) = ^elf)(fc)xa5(fc) . (B-27) 
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